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ABSTRACT 

In this paper I present the latest results on top quark physics from the 
D0 collaboration since the discovery of the top quark in March 1995. I 
summarize the discovery results, discuss progress since the discovery, and 
show how we can measure the top quark mass using three separate 
techniques. The measurements were made at the Fermilab Tevatron, a pp 
collider with V* = 1-8 TeV, using -50 pb ' data collected from 1992 to 
early 1995. 

1. Introduction 

Prior to the discovery of the top quark earUer this year, a surprising amount was 
known about its properties. A direct search for tt pair production by the CDF collaboration 
in 1992 using 4.1 pb ' of data^ led to a 95% confidence level (CL) lower limit on the top 
quark mass of 91 GeV. This limit was extended by the D0 collaboration using 13.5 pb"^ 
of data^ in 1994 to give a 95% CL lower limit of 131 GeV. In mid- 1994, CDF announced 
evidence for ti production from 19.3 pb ' data^ with 12 events found when only 5.7 ± 0.5 
were expected from known background processes. Since some of these events were 
double-Z7-tagged, this formed a 2.8 aeffect. In D0, we reanalyzed our 13.5 pb ' of data in 
Ught of our 131 GeV lower mass hmit, to search for higher mass top pairs''. We found 
9 events with a background of 3.8 ± 0.9 events, a 1.9 a effect consistent with no top 
production. These results were all in agreement with the indirect evidence for the existence 
of the top quark from fits to precision electroweak data from LEP, SLD and neutrino 
scattering experiments^. The fits at that time gave a mass window for the top quark of 
150 to 210 GeV, with the range determined principally by the choice of Higgs boson mass. 
So the stage was set for what was to come when D0 and CDF analyzed their much larger 
data sets from the Tevatron Run lb with -50 pb ' integrated luminosity. 

2. Top Quark Pair Production at the Tevatron 

The Tevatron Collider at Fermi National Accelerator Laboratory near Chicago in the 
United States is a proton-antiproton machine operating at a center of mass energy of 
= 1.8 TeV. At this energy, top quarks are produced predominantly in pairs, with -90% 
of the cross section coming from qq^ti interactions and the remaining -10% from 
gluon fusion gg^ti . In the analysis that follows, we assume that the top quark decays 
100% of the time as prescribed by the Standard Model (SM), i.e. t W^b and i W~b . 
This is a good approximation, since the CKM matrix element Vtb which links the top and 
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bottom quarks is very close to unity: 0.9988 < \Vtb\ < 0.9995 at the 90% CL^ assuming 
only three quark generations. We classify the various decay channels of the ti pairs by 
how the V7 bosons decay. When both Ws decay to an electron or muon (either directly or 
via a tau decay W — > tv-j:, t — > eVg or ixv^), then the ti channel is called dilepton 
production. When one W boson decays leptonically as above, and the other decays 
hadronically ( ud or cs and charge conjugate decays) then the mode is known as 
lepton+jets production. Finally, when both Ws decay hadronically, then this forms the 
alljets channel. Decay modes with a hadronically decaying tau lepton in the final state are 
not considered separately in this analysis, but form a small subset of the lepton+jets modes 
from Te, T^i and tt production, and a similar small subset of the alljets channel from 
T+jets and tt production. The branching fractions for the dilepton decay modes are -1.2% 
each for ee and fifx decays, -2.5% for the ep. mode, -15% for each of the lepton+jets 
modes (e+jets and ^u+jets), and -45% for the alljets channel, plus small contributions 
from the tau channels. 

The analysis strategy adopted by D0 for searching for ti events was to select events 
with isolated leptons, jets and missing transverse energy (evidence for one or more 
neutrinos in the event) and to determine the event selection efficiencies from data. We 
measured all the principal backgrounds and many of the minor ones directly from the data, 
and checked the background calculations using a looser event selection which allowed 
many more events to pass the cuts. Relative to the most recently published D0 results^ the 
search which led to the discovery had a data set approximately four times as large, the new 
analysis reduced the backgrounds by a factor of between three and six (depending on the 
decay channel) and the signal acceptance was kept between 70 and 80% of the previous 
value. 

3. The D0 Detector and Particle Identification 

The D0 detector is described in detail elsewhere^. For this analysis, the critical 
elements included a nonmagnetic central tracking system for tracking electrons, muons and 
jets and rejecting photons, a calorimeter for measuring the transverse energy of leptons and 
jets and the total missing transverse energy for inferring the passage of a neutrino, and an 
outer muon spectrometer for measuring the momentum of the muons. Particle 
identification algorithms are also described in detail elsewhere*^. Since there is no central 
magnetic field, D0 cannot distinguish between electrons and positrons and they will be 
known generically as "electrons" here. An electron is identified by the presence of an 
electromagnetic cluster of energy in the calorimeter with a track pointing towards it from 
the primary vertex. The energy cluster must be isolated from other clusters, and in the 
central pseudorapidity region of the detector. It must not be identified as a pion by the 
transition radiation detector. A muon is identified by a central track which connects the 
primary vertex to a track passing through the muon spectrometer, via a track through the 
calorimeter along which is deposited the energy from approximately one minimum 
ionizing particle. If the muon is considered to have come from a W boson decay, then it 
must be isolated from any jets in the calorimeter and from any other leptons in the event. If 
it is considered to have come from the semileptonic decay of a hadron in a jet (from the 
b quark from the top decay) then it must be close to or within the radius of a jet. It is then 
known as a tagging lepton, or "^u tag" and is used to identify a jet as a "Z? jet". 
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4. Data Search 

(a) Dilepton Channels 

The principal advantage of search for pp^ttX^ee, jifji, efx+X is that these dilepton 
channels are very clean, i.e. there are almost no background processes that mimic the 
signals. The disadvantage is that the branching fractions are very low. The cross section for 
ti dilepton events is only 3.5 x 10^^ of the total pp cross section, for a top quark mass of 
180 GeV. The main backgrounds are dilepton decays of Z bosons with accompanying 
initial state radiation to fake the b jets from top, Drell Yan production, W^W +jets, WZ, 
multijets with bb or cc decaying semileptonically, and instrumental fakes where a jet is 
misidentified as an electron and there is significant mismeasured missing transverse 
energy. Table 1 shows the cuts applied in the dilepton search for each channel. is the 
scalar sum of the transverse energies of the jets and the electron (if present) in the event. To 
study the backgrounds more carefully, the Hj cut is removed. The cuts are then known as a 
loose selection. Table 2 shows the results of the search for top pair production in the 
dilepton channels. 



Cut 


ep 


ee 


pp 


Lepton pj. [GeV] 


> 15/ 12 


> 20 / 20 


> 15/ 15 


Missing £y [GeV] 


>20 


>25 


>0 


Number of Jets 


>2 


>2 


>2 


Jet^r [GeV] 


> 15 


> 15 


> 15 


[GeV] 


> 120 


> 120 


> 100 


not aZ 






Z kinematic fit 



Table 1 Dilepton event selection cuts. 



Channel 


ep 


ee 


pp 


Total 


\Ldt [pb-^] 


47.9 ± 5.7 


55.7 ±6.7 


44.2 ± 5.3 




MC tt Signal 


0.34 ± 0.04 


0.25 ±0.05 


0.11 ±0.02 


0.70 ± 0.07 


Background 


0.12 ±0.03 


0.28 ±0.14 


0.25 ± 0.04 


0.65 ±0.15 


Data 


2 
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Table 2 Dilepton search results. The tt MC calculation is for a top quark of mass 
200 GeV, using the ISAJET generator. 



From approximately 50 pb~' of data, we found two ep events and one pp event, 
with an expected background of 0.65 ± 0.15 events. The ISAJET^ Monte Carlo signal 
prediction for the number of events expected is for a top quark mass of 200 GeV. Figure 1 
shows the distribution in of the three dilepton channels of ti production, and the 
respective backgrounds. It can be seen that the variable provides good discriminating 
power between signal and background. Hj may be considered to be a measure of the 
temperature of the event. 
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(b) Lepton+Jets Channels 

The main advantage of using the lepton+jets channels to search for tt production is 
that the branching fraction is significantly higher than for the dilepton channels. This 
advantage is counterweighed by the much higher backgrounds. The backgrounds to these 
channels are principally W+jets events and multijets where a jet fakes an electron and there 
is significant mismeasured missing transverse energy. 

We have used two independent analysis techniques to separate — > e + jets and 
H + jets from the background processes. We split the data set into two nonoverlapping 
sets, one set containing a tagging muon which was used to identify a jet, and the other set 
without such a muon, to which we applied a topological selection procedure. The variables 
used for event shape selection included Hj^ here defined as the scalar sum of the transverse 
energy of all the jets, and aplanarity^" A, a measure of how planar the event is. The tagging 
muon had to be nonisolated, i.e. close to or inside a jet, with the opposite definition to that 
used to label a muon as isolated (one from a W boson decay). Since requiring a tagging 
muon is quite a stringent demand, we relax the number of jets required from 4 to 3 to 
regain acceptance, remove the aplanarity cut completely and loosen up the Hj. cut. Table 3 
shows the event selection cuts used for lepton+jets events. The loose cuts involved 
removing the //^ cut and relaxing the aplanarity cut from 0.05 to 0.03 in the nontag event 
sample. Table 4 shows the results of these cuts on the search for tt pair production in the 
lepton+jets channels. 



Cut 


e+jets 


//+jets 


e+jets/ lii 


//+jets/ lii 


Lepton pj. [GeV] 


>20 


> 15 


>20 


> 15 


Missing E'^ [GeV] 


>25 


>20 


>25 


>20 


Number of Jets 


>4 


>4 


> 3 


>3 


Jet^r [GeV] 


> 15 


> 15 


>20 


>20 


Aplanarity (A) 


>0.05 


>0.05 


>0 


>0 


[GeV] 


>200 


>200 


> 140 


> 140 


Tagging [GeV] 


none 


none 


>4 


>4 



Table 3 Lepton+jets event selection cuts. 



Channel 


e+jets 


;U+jets 


e+jets/ in 


//+jets/ in 


Total 


jUt [pb-^] 
MC tt Signal 
Background 
Data 


47.9 ±5.7 
1.84 ±0.31 
1.22 ±0.42 
5 


44.2 ± 5.3 
0.95 ± 0.24 
0.71 ±0.28 
3 


47.9 ±5.7 
0.81 ±0.16 
0.85 ±0.14 
3 


44.2 ± 5.3 
0.41 ±0.10 
0.36 ±0.08 
3 


4.01 ± 0.44 
3.14 ±0.53 
14 



Table 4 Lepton+jets search results. The tt MC calculation is for a top quark of 
mass 200 GeV. 
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Background subtraction for the lepton+jets channels was done in the following 
manner. Two independent methods were used. The first invoked Berends Scaling^\ which 
notes that the ratio of the number of events with at least n jets to the number of events with 
at least n-l jets should be a constant for any number of jets. Therefore, the numbers of 
events with 1 or more jets and with 2 or more jets (after subtraction of the multijet 
background where there is no real W boson) are used to predict the numbers of events with 
3 or more jets and with 4 or more jets. The multijet background is subtracted first, because 
although it too should obey Berends scaling, the ratio between events with different 
numbers of jets is slightly different than that in the W+jets background. This technique was 
cross-checked with an event density fit in the {A-H^) plane and the results were found to be 
consistent with each other. After all cuts were applied to the /+jets events with no ^ tag, 
there remained 8 candidate events with a background of 1.9 ± 0.5 events. 

The idea behind looking for a muon in or near a jet to label it as a Z? jet is as follows. 
Every ti event has two & jets from the direct decays of the top quark and i antiquark. The 
b quark can decay semileptonically via b — > c//v^ to produce a final state muon and also 
the daughter charm quark can decay via c sjxv^ to give a second muon. Either or both 
of these oppositely charged muons may be detected and used to label the b jet. Using the 
branching fractions for these semileptonic decays, one can see that -44% of all ti events 
will have at least one muon in them associated with a b jet. D0's probability for 
reconstructing such a muon in or near a jet is -45%, leading to an overall result that -20% 
of aU ti events will have an observed muon tag. The corresponding probability to tag a 
random background jet (the mistag probability) is -0.4%. 

In the analysis, the tagging probability was measured as a function of the number of 
jets in the event, the transverse energy of those jets and the amount of missing transverse 
energy. This probability was used to calculate the number of expected background events 
in the /+jets/ ix channels. Backgrounds in these channels come mainly from V7+jets events 
and multijet events where at least one jet is from a & or c quark which has decayed 
semileptonically to give a muon. The probability for a particle in a jet to punchthrough into 
the muon spectrometer to fake a muon in the jet is negligible. The measured probability for 
tagging a jet with a muon was convoluted with the spectrum of VK+jets events where there 
is no tag to give the total predicted background from W+jets and multijet events with a 
tag. This technique was cross-checked using the dijet and y+jets data samples and 
found to give consistent results. After all the cuts had been applied to the Z+jets/^u events, 
there were 6 candidate events left, with a background of 1.2 ± 0.2 events. Figure 2 shows 
the number of events versus the inclusive jet multiplicity for the data after subtraction of the 
multijet events (i.e. nominally just W+jets events), for the W+jets events which have no 
^ tag convoluted with the measured tag rate, and the W+jets events with an appropriate 
amount of Monte Carlo ti production added. It can be seen that the data show a significant 
excess over the untagged W+jets events convoluted with the tag probability, which is 
compensated for by the addition of the Monte Carlo ti events. Also, the excess between 
data and Standard Model expectation without ti production increases as the number of jets 
goes up, in line with expectations for the high multiplicity top events. 
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(c) Data Search Summary 

The results of the D0 search for tt production in the seven separate decay channels 
are shown in Table 5 for both the standard set of cuts and for the loose set. The loose cuts 
were designed to let more background events into the sample to enable us to check our 
background calculations. It can be seen that the calculated number of background events 
does indeed increase in a similar manner to the number of events which pass the loose 
selection cuts. We calculate the probability that the background fluctuates up to give at least 
the number of candidate events seen, and find that it is only 2 x 10"^ (corresponding to 
4.6 a if the errors on the backgrounds are Gaussian). Also shown in the table are 
measurements of the tt production cross section from each of the sets of cuts. The two 
values are consistent with each other. 

Additional properties of the candidate events which indicate that the excess events are 
due to top quark pair production include the distribution of the 17 candidates between the 
7 channels. The distribution is consistent with the top hypothesis at the 53% CL, which is 
very high. We calculated the cross section from each type of decay channel with a top 
quark mass of 200 GeV and found consistent results between each channel, which may be 
interpreted as an additional piece of evidence that the excess of events seen are due to ti 
pair production. These results are shown in Table 6, together with the cross section values 
calculated as a function of the top quark mass, using the results from the standard event 
selection cuts. 





Standard Cuts 


Loose Cuts 


Dilepton 


3 


4 


Lepton+jets (shape) 


8 


23 


Lepton+jets (b tag) 


6 


6 


Total # of Candidates 


17 


33 


MC ti Signal 


4.7 ± 0.7 


6.3 ±0.9 


Background 


3.8 ±0.6 


20.6 ±3.2 


Probability 


2 X lO '^ (4.6 a) 


0.023 (2.0 a) 


aiti^X) [pb] 


6.3 ±2.2 


4.5 ±2.5 



Table 5 Summary of search for ti production in all seven decay channels. The tt 
MC calculation is for a top quark of mass 200 GeV, as is the measured 
cross section. 



Figure 3 shows the results of the D0 measurement of the top quark cross section as 
a function of the top quark mass, with a resummed next-to-leading-order (NLO) theoretical 
calculation of the cross section shown superposed for comparison'^. 

The results from the D0 discovery of ti pair production were announced on March 
2nd 1995 at Fermilab and appeared in Physical Review Letters on 3rd April 1995'^ A 
similar publication from the CDF collaboration appeared simultaneously''*. They measured 
the cross section as 6.8^|;^ pb for a top quark mass of 176 ± 8 ± 10 GeV. 
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tt Decay Channel Cross Section 
(m, = 200 GeV) [pb] 



Top Quark Mass Cross Section 
[GeV] [pb] 



Dilepton 

Lepton+jets (shape) 
Lepton+jets (b tag) 
All channels 



7.6 ±5.8 
4.9 ±2.5 
8.9 ±4.8 
6.3 ±2.2 



140 
160 
180 
200 
220 
240 



16.2 ±5.6 
10.8 ±3.7 

8.2 ±2.9 

6.3 ±2.2 
5.1 ±1.7 
4.3 ±1.5 



Table 6 // production cross section for each of the decay channels when m, = 
200 GeV, and as a function of top quark mass from all channels combined. 



5. What's New Since the Discovery of Top? 

(a) Overview of Work in Progress 

Since the discovery of the top quark by D0 in March 1995, D0 has approximately 
doubled its data set, from -50 pb"^ to -100 pb"'. Much work is going into analyzing this 
additional data and new results will be available shortly, but are beyond the scope of this 
paper. We have released preliminary results of a search for tt production in the alljets 
channel which will be discussed briefly here. We have also made a first preliminary 
measurement of the top quark mass using the dilepton events, which will be presented in a 
later section of this paper. New search techniques are in progress, including the use of 
neural networks and probability density estimators. And finally, we have begun a search 
for single top quark production'^. 

(b) Search in the Alljets Channel 

In the alljets channel, both W bosons from top decay hadronically, leading to a final 

state with at least six jets (bbqq'qq') where the q quarks are light. The advantage of 
searching in this channel is that it contains -44% of the total branching fraction. The 
disadvantage is that there is almost overwhelming background from multijet events with at 
least six jets for which there is no Monte Carlo simulation. The cross section for multijet 
events with >6 jets is -1000 x the ti alljets cross section. Event selection proceeded 
using topological techniques, coupled with b tagging. The goal of the search was to obtain a 
signal-to-background ratio of 1:1, with a few signal events remaining. 

To make the search, we required at least six jets in an event within the pseudorapidity 
region |77| < 2 and with Ej.> 15 GeV, but we relaxed the transverse energy requirement on 
the sixth jet down to £'j.> 10 GeV to maintain good efficiency. We then performed a grid 
search on several topological selection parameters such as aplanarity, centrality, the average 
number of jets in an event and the transverse scalar energy of all the jets except the two 
most energetic ones ( H^). 
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Preliminary results from the alljets search are shown in Table 7 for two sets of cuts, 
the standard set designed to maximize the signal-to-background ratio, and a looser set 
designed to allow in some background so that it can be compared with the background 
predictions. The next step in this analysis is to calculate a cross section using the results 
from the search when two b tags are applied. When only one tag is found, then the ti cross 
section is measured to be 2.2 ± 4.7 pb using the standard cuts and 5.5 ± 4.8 pb from the 
looser cut search, for a top quark of mass 200 GeV. These values are consistent with each 
other, with the theoretical expectation and with the results from the searches in the dilepton 
and in the lepton+jets channels. 





Standard Cuts 


Loose Cuts 


no b tag 






MC ti Signal 


10.4 


43.4 


Data 


102 


1483 


1 b tag 






MC ti Signal 


2.3 


9.5 


Data 


4 


50 


Background 


3.2 ±0.3 


41.7 ±1.4 


Cross Section [pb] 


2.2 ± 4.7 


5.5 ±4.8 


2 b tags 






MC ti Signal 


0.7 


1.4 


Data 


1 


3 


Background 


0.2 


1.2 



Table 7 Preliminary results from the search for tt production in the alljets channel. 

Predicted yields for Monte Carlo ti signal are for a top quark of mass 
200 GeV. 



6. Measuring the Top Quark Mass 

(a) Introduction 

We have used our candidate top events to measure the top quark mass using the 
decay kinematics. We have developed three techniques which I present here. The first uses 
the lepton+jets events from the time of the discovery and makes a fit using kinematic 
constraints. The second method uses the mass of the W boson reconstructed from the two 
jets into which it decayed as a calibration point for reconstructing the top quark mass from 
its decay products. This method is still in progress and has not yet lead to a new result on 
the top quark mass. The third technique uses the dilepton candidate events from a larger 
data set than that used for the discovery, and makes a fit to extract the top mass. This 
analysis leads to a new preliminary value for the top quark mass. 
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(b) Mass Fitting in Lepton+Jets Events 

In lepton+jets events, the process is pp — > ttX — > (blv) + {bqq') + X, so there are at 
least four jets in the final state, plus the lepton and some missing transverse energy from 
the passage of the neutrino out of the detector. The mass fitting technique is to make a two 
constraint (2C) kinematic fit, by smearing the energies of the jets and lepton (using 
measured energy resolutions) to force reconstruction of the two W bosons with 
^iv = ^qq' = f^w and then to force the top quark and t antiquark to have the same mass 
as each other nit = mi. Sixteen variables are used in the fit: missing transverse energy, and 
[E, !],(!)) for each of the five objects in the event (b,b,q,q',l). There is one unmeasured 
variable, the z component (along the beamUne) of the neutrino momentum. It is not 
possible to unambiguously reconstruct p^, since the underlying event from the proton- 
antiproton interaction deposits much energy in z which is not measured. For this fit, only 
the four highest jets with jryj < 2.5 are used. Attempting to combine additional jets with 
these highest energy ones, under the assumption that the extra jets might be final state 
radiation worsens the mass resolution. This is because although some of these jets are 
indeed final state radiation, they are sometimes recombined with the wrong jet, and also 
there is a significant amount of very high energy initial state radiation which should not be 
added in with the decay product jets from the tops. The jets were reconstructed using a 
narrow cone radius of 0.3, where the cone radius is defined as R- ^Arj'^ + . 
Corrections were applied to the jet energy for out of cone gluon radiation. The narrow cone 
was chosen to optimize the correspondence between reconstructed jet and initial state 
parton from Monte Carlo studies. 

Figure 5 shows distributions of the fitted top quark mass with (a) exact parton 
energies from the Monte Carlo, and no initial or final state radiation, (b) a full GEANT^* 
simulation of the D0 detector applied to the final state partons, thus smearing the energies, 
and occasionally merging or splitting the jets, (c) the additional detrimental effects of 
including final state radiation in the simulation and (d) the full simulation with both initial 
and final state radiation as well as energy smearing from a full detector simulation. In each 
plot, the shaded histogram shows the mass distribution when the correct jets are combined 
to reconstruct the W, t and i , and the black line histogram shows the mass distribution 
when the up to three combinations of jets with fit < 7 are chosen. It can be seen from 
Fig. 5 (a) that the intrinsic resolution on the top quark mass is rather good, at about 10 
GeV, and the mean fitted mass is close to the top mass used in the simulation. Choosing 
the wrong combination of jets broadens this peak slightly. In plot 5 (b), the resolution of 
the D0 calorimeter widens the mass peak more, and choosing the wrong jet combinations 

significantly broadens the peak and shifts it to a lower mean mass. In 5 (c) the addition of 
final state radiation to the simulation shifts the peak to an even lower mass as energy is lost 
from the reconstruction and sometimes a gluon jet is more energetic than a quark jet from 
the decay, and so is mistakenly used in the fit. Finally, in 5 (d) the effects of initial state 
radiation are shown. Using an initial state gluon jet instead of one from a top decay can 
increase the reconstructed top quark mass, and so a high mass tail is added to the mass 
distribution. 

The reconstructed or fitted mass for each top candidate event has to be converted into 
an actual top quark mass, by making a correction for the tendency for the mass peak to be 
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shifted by initial and final state gluon radiation. By measuring the mean fitted top mass as a 
function of the input Monte Carlo top mass, one can make this correction. The relation 
used to make this correction is shown in Fig. 6. 

The results of making kinematic mass fits to the lepton+jet event candidates is as 
follows. Of the 14 events which pass the standard cuts, 1 1 have at least four jets in the 
central pseudorapidity region and are successfully fitted. It is calculated that of these events, 
2.1 ± 0.4 are actually from background sources. The true top quark mass using these 
events is found to be 199 GeV. Of the 29 events which pass the loose cuts, 24 have at least 
four central jets and are successfully fitted. There are 11.6 ± 2.2 background events in this 
sample. The true top quark mass is found to be 199 GeV, in excellent agreement with the 
result from the standard cuts analysis. The statistical and systematic errors are smaller 
using the loose cuts, and the final result is therefore taken from that analysis: 

mt = 199 til (Stat) (syst) = 199 Ij^ GeV 

The systematic error is dominated by the 10% contribution from the uncertainty in the jet 
energy scale. 

The upper two plots of Fig. 7 show in the shaded histograms the distribution of fitted 
masses of the candidate lepton+jets events for both the standard and loose event selection 
cuts. The outlined bins show which events have a muon tag (4 events which pass both the 
loose and standard cuts). The jet closest to the tagging muon was forced to be labeled as a Z? 
jet when combining jets in the mass fitting. The dotted curve shows the best fit to the 
distribution for the + jets candidates and the dashed curve shows the best fit to the 
background. The solid curve is the sum of the other two curves and shows the best fit to 
the data. The position of the peak of the best fit curve to the signal does not change 
significantly if the number of background events is held fixed or is allowed to vary as a 
free parameter in the fit. 

The lower two plots in Fig. 7 show the negative log likelihood of the fit to the signal 
as a function of the true top quark mass. The fits to the signal shown in the previous two 
plots were repeated over the full range of input top quark mass values, and the quality of 
the fit evaluated at each point (in 10 GeV intervals). The fitted masses were than converted 
to true masses using the relations shown in Fig. 6 and the fit quality plotted versus this true 
mass. The minima of the parabola-like curves shown here correspond to the best 
measurements of the top quark mass, both at 199 GeV. 

(b) An Aside - Properties of the ti Events 

Now that we have reconstructed the tt pairs from their decay products, we can study 
properties of the tt system such as the invariant mass of the pair and the mean transverse 
momentum of the top quarks. There are many theoretical extensions to the Standard Model 
which affect these distributions and so they will become of considerable interest in the 
future. At the moment there are too few events to make any precise differentiation between 
the various models. The upper two plots in Fig. 8 show the invariant mass distribution of 
the top candidates for the standard and loose cuts. The data is shown in the black crosses. 
The darker gray shaded histogram shows the predicted signal and the lighter gray 
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histogram shows the predicted signal plus background. The data are consistent with the 
Standard Model prediction within errors. 

The lower two plots in Fig. 8 show the mean transverse momentum of the 
reconstructed tops in these events for both standard and loose cuts. As before, data is 
represented by the crosses and Monte Carlo calculations of the signal and 
signal+background by the shaded histograms. No anomalies are seen. 

(c) (m-m^) Analysis 

The motivation for developing the following technique for extracting the mass of the 
top quark from the candidate events is to see if we can find a mass peak for the W boson 
from the two jets it decayed into. We will then be able to calibrate the top quark mass 
against the known boson mass. This cannot be done using the 2C kinematic mass fitting 
method discussed previously, since explicit W mass constraints were applied. (The jet 
energies were adjusted until they reconstructed to give a perfect 80 GeV W boson.) The 
method discussed here will not constrain the Wmass. 

As before, only the four highest transverse energy jets are used in the reconstruction, 
this time using a cone of radius 0.5 to identify the jets. The wider cone was chosen so as 
not to be biased by the jet energy corrections, which had been tuned to reproduce the 
V7 boson mass. For each of the four permutations of how to combine the jets (two 
permutations if there was a b tag), we defined two top masses mf^ and m!^^ and a dijet 
mass m^. We assigned a weight to each permutation based on the difference of top masses 
between the leptonic and hadronic sides of the decay,l m^t'^-mf'"^\. We then plotted {nit) 
versus for each weighted permutation of every event. We did not adjust the dijet 
energies and the dijet pairs were not selected for consistency with the W mass. The results 
of this procedure are shown in Fig. 9. The upper two plots show the (m-m^) plane for (a) 
200 GeV mass top HERWIG" Monte Carlo and (b) background W+jets VECBOS" MC 
events and multijets data. A clear peak at (m, = 190 GeV, = 11 GeV) as expected for 
the top MC shows that the method is successful at reconstructing both the top quark and 
the W boson. The nine highest bins in this region are lighter shaded in the figure to 
distinguish them. The peak in this first plot is projected onto the m, and axes in the 
lower two plots. No such peak is seen in the distribution from the backgrounds in the 
upper right-hand Lego plot, where the jets are mainly from lower energy gluons. 

The data and best fit combination of backgrounds and ti signal with m, = 200 GeV 
are shown in the (m-m^) plane in Fig. 10. The same nine bins are lightly shaded as in 
Fig. 9 where the ti MC signal peaked, to guide the eye. One can clearly see that the data 
selected with the loose cuts is a mixture of the tt signal and the backgrounds, with peaks in 
both regions of the plane. 

So finally, we take the event distributions shown in Fig. 10 and project them onto the 
two axes to see if there is evidence for the decay of a V7 boson as indicated by a peak in the 
dijet mass distribution. We also look for a peak in the region of a heavy top quark on the 
other axis. These projections are shown in Fig. 1 1. On the upper plot, a cut has been made 
on events which have a dijet mass below 58 GeV, in order to remove background. A peak 
is seen in the reconstructed top mass distribution consistent with the top mass found from 
kinematic fitting. In the lower plot we make a cut on events with a reconstructed top quark 
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mass less than 150 GeV to reject background. We now see a clean peak in the data above 
the predicted background, centered just below the boson mass as expected. 

This analysis technique has not yet been used to measure the top quark mass from 
the reconstructed dijet Wmass. However, we can say that a W — > jetjet peak is seen in the 
data, which is expected if the heavy particles decaying are top quarks. The probability of the 
background fluctuating up to give the observed peak in the (m-m^) plane is only 1.3%. 

(d) Dilepton Mass Analysis 

D0 has made a preliminary measurement of the top quark mass using the dilepton 
events. This technique is based on one proposed by Kondo'^ and later by Dalitz and 
Goldstein". The method has been extended by D0 to use the missing transverse energy, 
and was published applied to one ti e/jX event from the first collider run^. For the 
analysis presented here, a data sample of -72 pb ' has been used. In this sample, 5 dilepton 
candidates have been found (one ee, two e/j. and one /j-jj. events) with an estimated 
background of one event. To make the mass fit, we assume the two b jets are the two 
highest transverse energy jets in the event. Monte Carlo ensembles of 5 event experiments 
are run to obtain the statistical error of the measurement. The systematic error as before 
comes mainly from the jet energy scale, and also from the choice of MC generator. The 
upper plot of Fig. 12 shows the best fit resolution functions for the data (solid curve) and 
the background (dashed curve), as well as the fitted mass for the five candidate events. The 
fit is repeated for many possible top quark masses and the log likelihood of the fits is 
shown in the lower plot. The maximum log likelihood occurs for a true top mass of 
145 GeV, where the best fitted mass of 162 GeV has been corrected to give the true mass 
in a similar manner to that shown in Fig. 6. The preliminary new result from mass fitting 
the five dilepton event candidates is: 

m, » 145 ± 25 (Stat) ± 20 (syst) = 145 ± 32 GeV 

If one uses only the two e/i events, which have the smallest chance of being background, 
then the result is almost exactly the same as when all five dilepton candidates are used. 

7. The Future of Top Physics at D0 

Top physics is being actively pursued at D0 in many directions. We are extending 
the analyses from the 50 pb ' of data presented here to use the -100 pb ' which we 
accumulated by July 1995. Results from this large data set should be forthcoming in 
January 1996. In addition to more data, we are also expanding our analysis techniques to 
gain higher efficiency for event selection and better background rejection. We are 
improving our lepton identification methods, and extending the b tagging method to 
encompass electrons as well as muons. Much work is also going into improving the mass 
measurement techniques. A search for electroweak single top production has started. 
Finally, we are actively preparing for the next run in 1999, which should bring us 2 fb"^ of 
data, together with a significantly upgraded detector. We will have a central 2 T magnetic 
field, and a silicon tracker vertex detector which will enable us to do detached secondary 
vertex tagging of B hadron decays in jets, as well as the semileptonic tagging we currently 
use. The prospects for top physics in the next run are very exciting. 
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8. Summary of Top Physics at D0 

The D0 collaboration discovered the top quark at the Fermilab Tevatron collider in 
March 1995. In a -50 pb~' data set, we found 17 candidate ti events with a background of 
3.8 ± 0.6 events using our standard event selection cuts, and 33 events with a background 
of 20.6 ± 3.2 events using a loose set of cuts. We measured the top quark mass using two 
separate techniques. From the lepton+jets events, we found nit = 199^|^ GeV, and from 
the dilepton events we made a preliminary measurement of ~ 145 ± 32 GeV. At the top 
mass of 199 GeV, we measured the pp — > tiX production cross section to be 6.4 ± 2.2 pb. 
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Fig. 1 Dilepton Hj^ distributions. The dark gray line at high Hj shows the Monte 
Carlo tt distribution (m, = 180 GeV), the light gray line shows the 
background. The arrows indicate the Hj^ cut value for each channel. The 
black rectangles show the candidate events and the hatched rectangle shows 
an ee event which fails only the Hj^ cut. 
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Fig. 2 Lepton+jets/muon tag events versus the inclusive jet multiplicity, before the 
Hj cut was applied. The pale gray points with large error bars show the data 

(with multijets subtracted). The mid-gray points show W+jets data where 
there is no tagging muon found, convoluted with the measured tagging rate. 
The dark gray points show how the W+jets/no tag x tag rate data become 
consistent with the total tagged data when ti Monte Carlo events are added. 
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Fig. 4 Distribution of events with at least six central jets in the vanauic nj 
(defined in text). The dark solid line shows the Run lb data (-25 pb '), the 
dashed line shows the Monte Carlo ti distribution form, = 200 GeV and the 
dotted line shows the distribution form, = 160 GeV. At large Hj , the signal- 
to-background ratio starts to improve significantly. 
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Fig. 5 Distribution of fitted top quark mass from 180 GeV ti ^ e + jets events from 
an ISAJET simulation, with loose cuts applied. The shaded histograms show 
the distributions using the correct jet combinations. The black line 
histograms show up to three low solutions trying all possible jet 
combinations. 
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. 6 Mean fitted top quark mass versus true top quark mass. The solid black line is for the loose 
cuts and the dashed line is for the standard cuts. The standard cuts gain extra background 
rejection in the low mass region, but also throw out more low mass top events, thus biasing the 
distribution. The thin dotted line is to guide the eye for a one-to-one relation. The gray lines from 
mean fitted masses of 182 GeV (loose cuts) and 187 GeV (tight cuts) show how both 
measurements give a true top mass of 199 GeV. 




Fig. 7 The upper plots show for standard and loose cuts the fitted top quark mass 
for the data (shaded histogram) and the best fits to the signal (dotted curve), 
background (dashed curve) and total (solid curve). The outlined histogram 
bins are four b tagged events. The lower plots show the negative log 
likelihood of the true mass for the standard and loose cuts, with minima at 
199 GeV on both plots. The width of the parabola-like curve gives the 
statistical error of the measurement. 
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Fig. 8 Upper plots show the invariant mass of the tt system with data (crosses), 
Monte Carlo signal, 200 GeV ISAJET simulation (dark histogram) and MC 
signal+background (light histogram). Lower plots show the mean transverse 
momentum of the top quarks in these candidate events. 
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Fig. 9 Upper plots show the distribution of Monte Carlo signal (left plot) and 
background (right plot) in the (m-m„) plane. Lower plots show the 
projections of the top MC signal onto the m„ axis and the m, axis. 
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Fig. 10 Data with loose cuts (upper plot) and best fit backgrounds plus 200 GeV 
HERWIG MC tt signal (lower plot) in the {m-m„) plane. A clear peak 
can be seen in the data away from the backgrounds. 




Fig. 1 1 Projections of the tt ^ I + jets data onto the reconstructed top mass axis 
(upper plot) and onto the dijet mass axis (lower plot) of the (m-m„) 
plane, with cuts applied to reject background. The dark gray histograms at 
the front are background, the light gray histograms are data and the mid 
gray histograms at the rear are 200 GeV HERWIG tt MC + background. 




Fig. 12 The upper plot shows distributions of fitted mass for the data (solid line), 
backgrounds (dashed line) and candidate events (arrows). The lower plot 
shows the log likelihood distribution of the fits, with the peak at a true top 
quark mass of 145 GeV. 



